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Abstract 

Total cross-sections and logarithmic slopes of the elastic scattering cross- 
sections for different hadronic processes are calculated in the framework of 
the model of the stochastic vacuum. The relevant parameters of this model, 
a correlation length and the gluon condensate, are determined from scatter- 
ing data, and found to be in very good agreement with values coming from 
completely different sources of information. A parameter-free relation is given 
between total cross-sections and slope parameters, which is shown to be re- 
markably valid up to the highest energies for which data exist. 
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1. INTRODUCTION 

It is certainly a great challenge to establish a microscopic description of high-energy 
scattering in the framework of the field theory of strong interactions, i.e. QCD. Sophisticated 
treatments of perturbation theory! have lead to interesting results, which however are either 
qualitative or not able to explain the most striking phenomena. There is a rich amount of 
data on soft high-energy scattering , i.e. elastic scattering at high energies and momentum 
transfers smaller than the hadronic scale (~ 1 GeV). The newest and most precise data come 
from proton-antiproton scattering^, extending up to center of mass energies s/s = 1 , 800 
GeV. There are older data in the pp,Trp and Kp channels!'!, and for other hadronic channels, 
such as the T,p system!, they are still scarce. More soft-scattering data are expected for the 
near future from Fermilab and LHC. 

The energy dependence, in the full range of available data, is well described in the Regge 
picture! . The total cross-sectionsB increase with energy like s 0808 ^ leading to a hypercritical 
pomeron intercept. The variation of the slope of the elastic scattering cross-sections is also 
well described in the Regge picture with a slope of the pomeron trajectory a'(t) = 0.25 
GeV" 2 . 

The value of about 2/3 for the ratio of np to pp (or pp) total cross-sections, as well as 
certain factorization properties, are suggestive of an additive quark model!, in which the 
main features of high-energy scattering can be described through quark-quark scattering 
amplitudes. On the other hand, there is also a remarkable flavour dependence of the cross- 
sections, which decrease with the increasing number of strange quarks in the scattering 
channel. Such a feature is most naturally explained in models in which the cross-sections 
depend on the sizes of the hadrons, which is also indicated by the dependence of the slopes 
of the elastic differential cross-sections on the hadron sizes!. 

In this paper we evaluate elastic scattering amplitudes of hadrons in the framework 
of the model of the stochastic vacuum (MSV), originally developed in order to treat non- 
perturbative effects in low-energy hadron physics0'0. The model therefore deals with pa- 
rameters of non-perturbative QCD that play an essential role both in hadron spectroscopy 



and in high-energy scattering. Our treatment is based rather on loop-loop than on quark- 
quark scattering. An important consequence is that the high-energy cross-sections depend 
on the sizes of the hadrons, and this effect is due to the same mechanism that leads to 
confinement. 

The paper is organized as follows. In sec. 2 we give the theoretical foundations of 
our model. In subsec. 2.1 we shortly recapitulate essential aspects of the analysis of soft 
high-energy scattering and in subsec. 2.2 some important features of the model of the 
stochastic vacuum are explained. In sec. 3 we apply the MSV to soft high-energy scattering 
and in sec. 4 we evaluate the formalism derived in sec. 3 and find convenient numerical 
representations for the results. In sec. 5 we discuss the choice of input parameters and 
compare our theoretical results with experiment. In sec. 6 we present concluding remarks. 
In the appendices we discuss some more technical points. 

2. THEORETICAL BACKGROUND 
2.1 Soft high-energy scattering in non-perturbative QCD 

The first field theoretical approaches to soft high-energy scattering at small momentum 
transfer were attempted by Amati, Fubini and Stanghelinilll in the framework of the multi- 
peripheral model and by Gell-Mann, Goldberger and Low in a massive vector-exchange 
theoryEl. Since when it became clear that QCD is the fundamental theory of the strong 
interactions, many efforts were made to explain the relevant features of soft high-energy 
scattering, either by assuming genuine non-perturbative effectsB or by elaborate summation 
of perturbative contributionsS . 

The success of methods, applied primarily in hadron spectroscopy, taking into ac- 
count non-perturbative contributions as gluon-condensates0 stimulated Landshoff and 
NachtmannS to apply non-perturbative concepts also to soft high-energy scattering. In 
an Abelian model they related elastic high-energy scattering to the nontrivial structure of 
the QCD vacuum, requiring, besides the (static) gluon condensate, a finite correlation length 
for the slowly varying (non-perturbative) gluon fields in the vacuum. In their Abelian model, 
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this correlation leads effectively to a non-perturbative gluon propagator. The main conse- 
quences of the model were the correct spin structure (vector-like exchange) of soft elastic 
high-energy scattering and quark additivity. The latter holds only if the above mentioned 
correlation length of the gluon fields is small compared with the hadron radius. The energy 
dependence of the total cross section comes out only approximately correct, namely as a 
constant, instead of the experimentally observed slow rise like s - 0808 . This model has been 
applied successfully also to other channels^. 

The structure of non-perturbative contributions to high energy scattering was further 
investigated in a more general way by NachtmannEl, who reduced the non-perturbative 
parts in soft high-energy hadron scattering to quark-quark (anti-quark) scattering, justifying 
this reduction in certain kinematical regions through considerations in the femto-universe. 
Through the use of the eikonal method, the high-energy scale was separated and the (non- 
perturbative) part of the quark-quark scattering amplitude that is due to exchange of states 
with the vacuum quantum numbers could be reduced to an expression with the structure 

(q C '(pi) q D '(p'2)\q C \Pi) <?> 2 )> ^ 

U(p[) 7 M «(pi) U(p' 2 ) 7m M (^2) Jqq{q 2 ) Sc'C^D'D , (1) 

with small momentum transfer q = p[ — p\ = p' 2 — p 2 . The upper indices C, D denote the 
colours of the quarks, and u(p) is a Dirac spinor. 

The transition from the quark-quark scattering amplitude to the observable hadron- 
hadron scattering amplitude was achieved through current matrix elements occurring in 
deep inelastic scattering. 

In this paper we widely follow the general analysis of NachtmannEl, and evaluate quark- 
quark scattering amplitudes in a specific non-perturbative model, namely that of a stochastic 
vacuum with Gaussian fluctuations^ of the field strength. There is however a very specific 
difference: whereas the original treatment of Nachtmann is based on a reduction of hadron- 
hadron scattering to quark-quark scattering, our basic entities are scattering amplitudes 
for Wilson loops in Minkowski space-time. A definite advantage of our approach is the 



gauge invariance of these amplitudes in contradistinction to the quark-quark amplitudes. 
The loop amplitudes treated here can also be obtained in the framework of NachtmannEl, 
if one starts with hadron-hadron rather than with quark-quark scattering matrix elements 
(O. Nachtmann, private communication). 

Another important difference must be remarked. Quark additivity appears in a natural 
way in the Landshoff-Nachtmann mo del0 and in the extended framework of NachtmannEl. 
Here each quark interacts with the vacuum field, and we may consider that this interaction 
defines a region with the form of a tube around the quark path. Since the interaction with 
another quark is due to the correlations of the fields in two such tubes, the effective radius 
of a tube is actually determined by the correlation length of the vacuum field. The qq 
interactions occurs only in the regions where two such tubes overlap. If the separation of 
the quarks inside a hadron is large compared with the correlation length a, the interaction 
regions for the different pairs of quarks when two hadrons collide are indeed well separated 
from each other, and quark additivity holds. However, such an argument is dangerous in 
a non-perturbative treatment, since, for instance, it would not lead to the area law for the 
Wilson loop. Indeed, it is well known that even a short range correlation of the fields can 
lead to long range effects due to potentials, since the step from the fields to the potentials is 
essentially non-local. A good example of this phenomenon is the Bohm-Aharonov effectS, 
where the phase of an electron can be influenced by a magnetic field located far away. 
These long range effects might spoil quark additivity, and we show in sec. 2.2 that this is 
indeed the case with the model of the stochastic vacuum : the same effect which leads to 
confinement also leads to a violation of quark additivity. The physical reason for that is 
easy to understand, since not only the quarks but also the glue between them participate 
essentially in the scattering process. We return to this point in technical detail in sec. 2.2. 

2.2 The Model of the Stochastic Vacuum 

The model of the stochastic vacuum0'0 is based on the idea that the low frequency 
contributions in the functional integral can be taken into account by a simple stochastic 
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process with a converging cluster expansionE2l. This assumption leads, in a non-Abelian 
gauge theory, to linear confinement of static colour sources. 

Let us phrase this idea in a somewhat more formal language. Let V(f)(k) e~ s ^ be the 
functional measure of a quantum field theory where the fields to be integrated over are 
expressed in momentum space. If the measure is split into low and high frequency parts 

vm = n mk) n v~m , (2) 

the integration over the high frequencies can be accounted for perturbatively, in an asymp- 
totically free theory. All terms of higher than quadratic order in the fields occurring in the 
exponential of the action e~ s ^ are expanded in a power series and, since only quadratic 
terms are kept in the exponent, the remaining functional integrals are Gaussian. 

Little is known about the functional integration over the low frequencies. In the QCD 
sum rule approachEl, the contribution of this part is taken into account by power corrections 
proportional to specific non-perturbative vacuum expectation values (condensates). A model 
which goes further is that of the stochastic vacuum00 in its most restrictive form. Since we 
know that nature has managed to regularize the infrared problems of perturbative theory 
(after all, we do observe hadrons), we may assume that the integration measure of the low 
frequency fields may be approximated by a simple functional measure. The simplest ansatz 
is that of a Gaussian integration measure, which is specified by a correlator (corresponding 
to the propagator in perturbation theory). This correlator is (apart from its specific form) 
determined by two scales: the strength of the correlator and the correlation length. 

Generically, we may write 

/ J] V~${k)e- S <f>( Xl ) 0(x 2 ) = (<j>( Xl )<l>(x 2 )) A = G{ Xl - x 2 ) , (3) 
\k\<ii 

and obtain all other Green's functions for a Gaussian process by factorization. Note that 
we have liberally switched between the fields in coordinate space, <f)(x), and in momentum 
space, <p(k). Eq.(|3]) describes only the low internal frequencies part of the correlator, and thus 
G(xi — £2) is supposed to be regular for x\ — > x 2 . The singularities are due to perturbative 
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terms. As mentioned above, this simple model leads to confinement in a non-Abelian gauge 
theory, and moreover the heavy quark potential deduced from the correlation determining 
the Gaussian stochastic process agrees very well with phenomenological determinations^!!. 
Thus it is not unreasonable to apply the approximation by a Gaussian process to other 
non-perturbative phenomena, as soft high-energy scattering. 

First we discuss some basic properties of the two point correlator defining the non- 
perturbative Gaussian process (i.e. the approximation for the measure of functional inte- 
gration over the low frequency fields). As always occurs with approximations, a control of 
gauge invariance is essential. If that control is missing, one is never sure whether the output 
of the calculation is only a gauge artefact. So we do not deal with the correlator of the 
gauge potentials A?(x), but rather of the field-strengths F^ u , where the upper index F is 
the colour index of the adjoint representation. As usual, we introduce the Lie-algebra valued 
field quantities 

V 2 -l 7V C 2 -1 

A M (x) = £ r F , F^(x) = £ F% r F , (4) 

F=l F=l 

and the covariant derivative 

D " = - ujA ) ■ (5) 

In the expression above, Tp represents the N^ — l generators of the Lie algebra of the gauge 
group SU(N C ). For N c = 3 they are, in the fundamental representation, 1/2 times the Gell- 
Mann matrices \p. If not stated otherwise, we always use the fundamental representation. 
For some later discussions (see eq.(^) below), it is convenient to work with a general number 
of colours N c . In the first part of this section we work, as usual in the functional approach, 
in an Euclidean space-time continuum, and therefore only lower Lorentz indices are used. 

In a non-Abelian gauge theory the field strength tensor ¥^ v (x) does change under a local 
gauge transformation 

F^x) -> U(x) F^(x) V-\x) , (6) 
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where U(x) is a local element of the gauge group SU(N C ). 

In order to give a well defined meaning to a correlator, which is a bilocal object, we 
parallel-transport the colour-content of all fields to a single reference point w, i.e. we consider 
the parallel-transported field strength tensor 



F^ u (x;w) := 1 (x,w) F^(x) 



where <p(x, w) is a non-Abelian Schwinger string from point w to point x 



(7) 



(x, w) = P exp 



-ig j da (x — w)^A^{w + a(x — w)) 



(8) 



P denotes path ordering, which is necessary in order to give to the exponential a well defined 
meaning. In a non-perturbative way, this ordering is defined, for any operator O through 



P exp 



0(a) da 



Km TJexp 



iO( ) ■ Aa k 



(9) 



with < cti < 2 . . . < 1 , Aa k = a k+1 - a k . 

The field strength tensor in eq. (0) transforms with the gauge transformation at the fixed 
reference point w 



F^(x,w) -> U(w) F^(x;w) V~ L (w) . 



(10) 



The correlator (F flv (x, w) Fs a (y,w)) A , i.e. the vacuum expectation value with respect to 
the low frequencies , is gauge- covariant because of eq.fllUp; in general, it may depend on the 
two coordinate differences (x — w) and (y — w). We now make the crucial approximation 
that the correlator is independent of the reference point w, and thus only depends on the 
difference z = x — y. This approximation becomes exact for x — > y, and is nearly true for 
w fixed and large distances z = x — y. In this approximation, the most general form of the 
correlatoi0 is given by 



g 2 Fg(x,w) FZ(y,w) 



8 CD 1 

hg 2 ff) 



a AT2-1 12 
- V 5 vp ) ■ D(z 2 /a 2 ) 











D l (z 2 /a 2 ) 



(11) 



Here z = x — y, a is a characteristic correlation length, (g 2 FF) is the gluon condensate 

(g 2 FF) = (g 2 i£(0) F°(0)) A , (12) 



N c is the number of colours, C, D — 1, . . . , N 2 — 1 , and the factors in eq. (|TT|) are chosen in 
such a way that 

D(0) = A(0) = 1 . (13) 

The two possible tensor structures are arranged in order that the second term satisfies 
the homogeneous Maxwell equation, i.e. 

-^T £ ^(g 2 F°(x, w) F^ a {y, w)^ = for k = . (14) 

Hence, in an Abelian gauge theory without monopoles, where the homogeneous Maxwell 
equations must hold, only the second structure can occur, i.e. we must there have k = 0. 
However, in a non-Abelian theory there is no reason for k to be zero. 

With the form of eq.flTTl) for the correlator, one obtains00 the area law for a Wilson 
loop with the string tension p, given by 



K7V 



oo 



p = Tu {92FF)a J D ^ du2 - ( 15 ) 

Thus only the tensor structure proportional to D leads to confinement. This result has 
the very welcome consequence that only in non-Abelian gauge theories the model of the 
stochastic vacuum leads to confinement. It also teaches us that the correlator D is the part 
specific to non-Abelian gauge theories. 

The correlator in eq.(|Tl"|) has been calculated on the lattice^, and the results show 
unambiguously that k is different from zero, as predicted by the model of the stochastic 
vacuum. The ratio k/(1 — k) is rather large (about 3), so that D(z 2 /a 2 ) is the dominant 
contribution. We return to this point in more detail later. 

3. THE MODEL OF THE STOCHASTIC VACUUM IN HIGH-ENERGY 
SCATTERING 



The correlator in eg. (|TTD , specifying the Gaussian process that approximates non- 
perturbative effects of QCD, is the starting point for our evaluation of observables in soft 
high-energy scattering. In the analysis mentioned in sec. 2, Nachtmann0 evaluated the 
quark-quark scattering amplitude using the eikonal approximation for the interaction of the 
quarks with the gluon field. In a first step, we follow the same approach, and consider the 
scattering amplitude of a single quark in a given external colour potential A^. If the energy 
of the quark is very high and the background field has only a limited frequency range, the 
quark moves on an approximately straigth light-like line and the eikonal approximation can 
be applied. At the end of this section we recall the condition for the validity of the eikonal 
approximation. 

Along its path T, the quark picks up the eikonal phase (which is here a unitary N c x N c 
matrix) 

V = P exp[-ig J A^z) dz^\ . (16) 

Here A M is again the Lie- algebra valued potential and P denotes path ordering (see eq.([|)). 
The phase factor for an antiquark is obtained by complex conjugation. 

From the scattering amplitudes for single quarks in the background field, we obtain 
the non-perturbative quark-quark scattering amplitude by functional integration over the 
background field of the product of the two scattering amplitudes. More specifically, consider 
two quarks travelling along the light-like paths I\ and T 2 given by 

r x = (x°, 6/2, x 3 = x°) and T 2 = (x°, -6/2, x 3 = -x°) , (17) 

corresponding to quarks moving with velocity of light in opposite directions, with an im- 
pact vector 6 in the x 1 a; 2 -plane (referred to in the following as the transverse plane). Let 
Vi j2 (±6/2) be the phases picked up by the quarks along these paths 



^1,2 (±6/2) = Pexp 



-ig [ A^z) dz^ 



(18) 



Then the scattering amplitude for two quarks with momenta pi, p 2 and colour indices ci, c 2 
leading to two quarks of momenta p^, p^ and colours C3 C4 is given by0 
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T c 3Ci ;c lC2 (s,t) = U(p 3 ) 7 M 0l) W(P2) IMP*) V > ( 19 ) 

where 




Here ( ) A denotes functional integration over the background field; q is the momentum 
transfer (p\ — p 3 ) projected on the transverse plane. Of course the approximation makes 
sense only if |g| <C \p\. The quantity Z$ is the fermion wave-function renormalization 
constant in the eikonal approximation, given byEl 

z^A] = itr [y x (o)] = i-tr [y 2 (o)] . (21) 

The subtraction of the unit operator from the phase-matrices V is due to the transition from 
the S to the T-operator. 

In the limit of high energies we have helicity conservation 

u(p 3 ) 7Mpi) w(P4) 7mm(p2) s -^2s<5 Ai a3^A4A 2 , (22) 

where A« are the helicities of the quarks and s = (pi + P2) 2 - In the following we can thus 
ignore the spin degrees of freedom. 

The scattering amplitude in eq.(|l9|) is explicitly gauge dependent and the cautioning 
remarks made in the last section apply here. But we know that, in meson-meson scattering, 
for each quark there is an antiquark moving on a nearly parallel line. Furthermore, the 
meson must be a colour singlet state under local gauge transformations. To construct such 
a colourless state we have to parallel-transport the colour content from the quark to the 
antiquark (or vice versa) in the same way as discussed in sec. 2 for the field-strength tensor. 
Since this parallel-transport of the colours is made by a Schwinger string 4>(x g , Xg) (see 
eq.(@|)), we obtain for the meson a Wilson loop whose light-like sides are formed by the 
quark and antiquark paths, and front ends by the Schwinger strings (see fig. 3.1). The 
direction of the path of an antiquark is effectively the opposite of that of a quark, so that 
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the loop has a well defined internal direction. The resulting loop-loop amplitude is now 
specified, not only by the impact parameter, but also by the transverse extension vectors. 
We thus introduce the loop-loop scattering amplitude 



J(b,R u R 2 ) 



(itr Wi(0, i?i)) A (^tr W 2 (0, R 2 )) A } _1 



/— tr 

\N r . 



I ! 2 ( --,i?2 



1 



-tr 



(23) 



where b/2,Ri) is the Wilson loop 



Hi l-^Ri 



P exp 



-«o 9 Ajz) dz^ 



(24) 



The closed loop dS\ is a rectangle whose long sides are formed by the quark path Y\ = 
(x , 6/2 + Ri, x 3 = x Q ) and the antiquark path rf = (x , 6/2 — x 3 = x ) and whose front 
sides are formed by lines from (T, 6/2 + Ri,T) to (T, 6/2 — R\,T) for large positive and 
negative T (we will then take the limit T — > oo). 1^2(6/2, i? 2 ) is constructed analogously. 
The first factor in eq. (|23|) is the loop renormalization constant that replaces the quark field 
renormalization in eq. (f20|) . 

Our next aim is to perform the functional integration over A by applying the model 
of the stochastic vacuum discussed in the preceding section. Since the correlator is given 
in terms of the parallel-transported field tensor F IM/ (x,w), we have first to transform the 
line integrals J A^dz^ through integrals over the field tensor. This is done with the help of 
the non-Abelian Stokes-theoremB. By deforming the path as indicated in fig. 3.2, we can 
express the line integral § 9S A M (z) dz^, where the closed path dS goes from w to then 
to X{2) and back to w, into the surface integral f s F^ u (z,w) dE^fe), where dT,^ u (z) , with 
// < v, is the element of the surface S at point z. Here we have used that, for a sufficiently 
small contour, we have 



/ Apdz" = I V^dST + 0(S 2 ) 
J os J s 



(25) 



In this way we obtain 
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Pexp f -igA^z) dz» = P 5 exp / -igY^w) dW v {z) , (26) 

J dS J S 

where P5 now denotes surface ordering according to fig. 3.2. Since the reference point w in 
the correlator ( |TT] ) must be the same for both fields, we have to choose a common reference 
point for both traces in the product 



tr 



W 1 {- h -,R 1 )-l tr W 2 Q,R 2 )-1 



(27) 



We choose the point w in the most symmetric way and then the surface emerging from the 
loop dSi is formed by the sliding sides of a pyramid with the loop dSi as basis and the point 
with coordinates w as apex; the same holds for dS 2 (see fig. 3.3). 

Before its application to high-energy scattering, the model of the stochastic vacuum 
must be translated from Euclidean space-time, in which it is naturally formulated, to the 
Minkowski continuum. Unfortunately we cannot go the other way and continue eq.(f2~3|) 
to the Euclidean continuum, which would be the safe way from the point of view of the 
functional integration. However the Wilson loops occurring in eq.(p3|) have light-like sides 
which would shrink to a point if continued to a space time continuum with Euclidean metric. 
We think that this is a serious obstacle in all attempts to evaluate soft high-energy amplitudes 
numerically on a lattice. 

Since we cannot adapt the scattering amplitude to the Euclidean world, we have to 
proceed the other way and adapt the model of the stochastic vacuum to the Minkowski 
world. We are fully aware that this is by no means a trivial step and, pending a better 
analytical understanding of non-perturbative effects, we have to let the experiment decide 
on the justification (Similar problems occur when one applies instanton effects to high-energy 
scattering^!). Thus we must translate the correlation function in eq.([ll|) to the Minkowski 
world. This is obvious for the tensor structure, where we just substitute 5^ u by — g^ u , etc, 
but simple choices for the correlation functions like exp(— z 2 /a 2 ) or exp(— y/z^/a) cannot 
be analytically continued in a meaningful way. Therefore, we must look for correlation 
functions D(z 2 /a 2 ) and Di(z 2 /a 2 ) that fall off for negative z 2 values (corresponding to 
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Euclidean distances), and whose Fourier transforms exist in Minkowski metric, since these 
will enter essentially in the scattering amplitudes. An ansatz for the correlator that fulfills 
this requirement can be written in terms of the Fourier transforms 



6 CD 1 , r d 4 k 

{(flWSW ~ 9pa9up) KlD{k 



2\ 



where 



iD(k 2 ) = J d 4 zD(z 2 /a 2 ) e ikz/a (29) 

and 

iD^k 2 ) = J d 4 zD 1 (z 2 /a 2 ) e ikz/a . (30) 

After this choice is made, all functional integrations can be performed, in principle. 
The quantities W\p in eqs.(|23"D and (|4j) can be expressed as surface integrals, according 
to eq. (|26|) . The exponential being expanded, the expectation value can be calculated using 
eq.(|26|) and assuming factorization in a Gaussian process (see below). 

Before we enter into details, we make two remarks which facilitate further calculations. 
First we note that, since 

(tr W 1 ) A = (tr W 2 ) A = N c , (31) 

the functional integral over the surface of one pyramid alone vanishes. To see the formal 
reason for eq.([H]), we remark that in the evaluation of a single loop, say (tr Wi) A , only 
the expectation value (e+F lli (x ) w)e p + Fpk{y ) w)) A occurs, where e+ is the light-like vector 
(1, 0, 0, —1), and i, k (1 or 2) are indices of the transverse plane. These correlators are zero 
by virtue of the tensor structures given in eq . (p8|) (note that e+e +A1 = 0). Therefore only 
the unit term contributes, leading to eq . (|3lD . This is not in contradiction to the area law in 
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Euclidean space-time, since the area of loops with light-like sides would be zero in Euclidean 
metric. 

Due to this mechanism, the quantities 

Z lp = {^txW 1 (0,R i )) A (32) 



that enter in the first factor of the expression in eq. (|23[) are equal to one (no loop renor- 
malization). In the same way only mixed terms, from different pyramids, contribute for 
the expectation value in eq.(|23"D, since the correlation functions arising from the expansion 
of tr Wi(— 6/2, Ri) alone, or equivalently of tr ^2(6/2, R 2 ), contain only field projections 
etF^ or e^F^ respectively. 

We next expand the exponentials Wt in eq. (|23|) . Since in the expansion of the trace of 
the exponential at least two terms are necessary (tr = 0), and because of eq.(^TJ) the 
lowest order contribution to the loop-loop scattering amplitude is given by 

J(6,#i,# 2 ) = -H#) 4 (^[) 2 tr l T Ci T C 2 } tr [TDiTDj,] ■ 

2 



1^(f^ v1 { Xi ,w) F% U2 (x 2 ,w) F^ ai ( yi ,w) F°l 2 (y 2 ,w 



A 

+higher correlators . (33) 

We next apply the factorization hypothesis 

(F Cl F° 2 F Dl F D ' 2 )= (F Cl F° 2 )(F Dl F D ' 2 ) 

+ (F Cl F Dl ) (F c ' 2 F° 2 ) + (F Cl F° 2 ) (F° 2 F Dl ) , (34) 

where the arguments and the Lorentz indices of F Ci , F Di are the same as in eq. (|3"3|) . 

We have checked that the higher order terms are indeed small as compared to the leading 
term, and therefore we neglect them in the following. In this way the surface ordering 
indicated in eq.(^6|) becomes irrelevant. 

It is convenient to introduce the eikonal function \ 
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x(b,Ri,R 2 ) = {-igf I d£P»{x) I d^(y)(F^(x,w)Fg(y,w)) A . (35) 



' Si J S2 

Then the loop-loop amplitude J (6, R±, R 2 ) is given to the lowest order in the correlator by 



- - - 11 1 

J(b,R 1 ,R 2 ) = -—- 



X(b, Ri,R 2 ) 



(36) 



N 2 8 (N 2 - 1) 

We notice the presence of the colour suppression factor l/(iV| — 1) that always occurs in 
interactions between colourless objects. The eikonal function is determined by the geometry 



and by the correlator (28 



For the Fourier transform D(k 2 ) of the scalar correlation function D, that enters in 
eq . (|28D , we introduce an ansatz which fulfills the requirements made above, and write 

-6iA n k 2 1 



D(k 2 



{k 2 - l/\l)n \2n-6 



n > 4 



(37) 



with 



D{z 2 /a 2 ) = f D{k 2 ) e~ ikz/a 



d A k 



(38) 



These functions are discussed in Appendix 1. The constants A n and A ra are fixed by the 
normalization conditions 



POO 

D(0) = 1 and / du D(-u 2 ) = 1 

J 



(39) 



The second of these conditions allows the identification of a as a correlation length. 

The string tension for the correlator of the form (^) can be obtained^ from eq.(|l5|). 
With our choice for the correlator, we have§ 



p = K(g 2 FF)a 2 -(n-3) 



r(w-3) 
T(n - 5/2) 



(40) 



The scalar function D\ is completely independent from D, and may have different values 
for the parameters a and n. Lattice calculations^ show however that the forms of D and 
Di in the Euclidean region at large distances are similar, with D about 3 times larger 
than D\. We show in Appendix 2 that even with equal weigths for the two functions, the 
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contribution of D to high-energy soft scattering is by far the dominant one, so that we can 
safely concentrate on this function. 

Since analytic calculations are most easily done for the case n — > oo, previous calculations 
have been made with that choiceiHi, but for n = 4 there is good agreement with the 
form (exponentially decreasing at large distances) indicated by lattice calculations in the 
Euclidean regionil; we therefore adopt this form in the present work. 

Up to now we have only considered loop-loop scattering amplitudes. It would be highly 
desirable to have a formalism relating these fundamental field theoretical entities to observ- 
ables, corresponding to the operator product expansion for quark (and gluon) amplitudes. 
In the absence of such a formalism, we have to rely on a rather simple minded quark model. 

In a relativistic quark model, the distribution of the quarks is described by the transverse 
momentum k±_ and the fraction x of the longitudinal momentum carried by the quark. Since 
our amplitude is independent of the momentum of the quarks (as long as the energy is high 
enough to ensure light-like paths), we may neglect the dependence on x, and only consider 
the transverse dependence. This transverse dependence is given by the Fourier transform 
of the transverse wave-function, which determines the width 2|i?| of the Wilson loops. We 
thus obtain our hadron-hadron (here still meson-meson) scattering amplitude by smearing 
over the values of Ri and R 2 in eq.([3~3"l) with transverse wave-functions if>(R). This leads to 
the meson-meson scattering amplitude 

JMM>(b) = J d 2 Ri J d 2 R 2 J(b,R 1} R 2 ) \^ M (Ri)\ 2 \ipM>(R 2 )\ 2 ■ (41) 

The choice of the transverse wave-functions iPm{R) will be discussed in the next section. 

For the treatment of the baryons we restrict ourselves to N c = 3. We adopt two pictures: 
a genuine three-body configuration, and a diquark picture. In the latter the baryon is 
described exactly as a meson, where the diquark replaces the antiquark. In the three body 
picture the baryon is described as shown in fig. 3.4. There are three quark paths leading 
from X(i) to x'A), i — 1, 3. The coordinates Xh and x' h refer to the central point of the baryon. 
The paths from Xh to x^) and x' h to x',^ respectively must ensure that the baryon is a colour 
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singlet under local gauge transformations. This is done by parallel-transporting the colour 
from the quark positions xd) to Xh, and coupling the colours antisymmetrically in the form 

~^ £abc ( Paa'{x h ,X(i)) (f)bb'(x h ,X ( 2)) <Pcc'{x h , X(^)) , (42) 

where 0's are the Schwinger strings, eq.(§). An analogous factor occurs at the end (primed 
coordinates), so that the baryon is described by the product of paths 

^abc 4>ad[Fl\ 0f> e [r 2 ] C /[F 3 ] E de f , (43) 

where the path Tj leads from xh over x^) and x',^ to x' h . Let To be the path leading from 
x' h to Xh- Since cf)[T ] G SU(3), we have 

£def = Ba'b'c' 0da'[ r o] 4>eb'[^o] 4>fc'[^o] ■ (44) 

Inserting this expression into eq. (f43"D , we obtain that the baryon is represented by the product 
of Wilson loops (without traces) 

~ BabcW^idS,] W w [dS 2 ] W cc ,[dS 3 ] e a , Vd , (45) 
6 

where W / a a'[^S'i] = a <f[ri] 0da'|To] is the loop from Xh to xi, x[, x' h and back to Xh- Thus for 



a baryon the factor (l/JV c )tr [Wi(—b/2 } Ri) — 1] in eq. (|23[) has to be replaced by 



^ eabcBa'b'c' {w a , a (-&/2, Rj W Vh {-b/2, R 2 ) W dc {-b/2, R 3 ) - 8 aa , 5 W 5 CC ) . (46) 

Here Ri is the vector extending from the middle line T to the border of loop %. The impact 
parameter vector b is taken with respect to the middle line r . The factor 1/36 is due to 
colour normalization. We discuss baryon transverse wave-functions in the next section. 

We conclude this section with some considerations concerning the validity of the under- 
lying quark model, reviewing the analysis made by NachtmannS, who discusses limits on 
the energy values y/s in order to ensure the validity of his treatment. It is clear that the 
eikonal approximation can be justified only for a sufficiently high energy. A detailed analysis 
yields the condition y/s > 2r h~ 2 where h is a typical hadronic scale (h ~ 1 GeV^ 1 ) and 



To is the observation time of the scattering process in the femto universe. There is another 
scale Qo that indicates the separation between the perturbative and the non-perturbative 
effects. The correlation length and the wave-function depend on its value, since it can be 
viewed as a renormalization scale for non-perturbative quantities. The value of Qo should 
also be approximately equal to the hadronic scale h" 1 , and its relation with the energy and 
the observation time was found@ to be 

Ql « ^ ■ (47) 

From this relation Nachtmann also obtained a conservative upper bound for the energy, 
which results from the requirement that the observation time should be short enough so 
that the string between the scattering quarks or antiquarks is not broken if the two particles 
fly apart. The string breaks if the two quarks become so far apart that the potential energy 
of the string reaches the proto-hadron mass m c « 1.3 GeV . Then the breaking time r cr is 
determined by 

m c , „„s 

Tar = — , (48) 
P 

where p is the string tension p rs 0.18 GeV~ 2 . The maximal observation time is given by 
To = 2r cr . In the present paper we consider only elastic scattering of loops, where no string 
between quarks of different hadrons is formed, hence we see no compelling reason to apply 
the above limit to r , and there is no upper limit on the energy for the applicability of the 
model from this consideration. However, for defmiteness, we take in the present work as a 
reference energy the value yfs = 20 GeV, corresponding to the conservative upper bound 
indicated by NachtmannB, and discuss other energies separately. 



4. NUMERICAL EVALUATION AND PARAMETERIZATION OF THE RE- 
SULTS 

We now introduce the notation R(I, L), where the first index (1=1,2) specifies the hadron, 
and the second specifies the particular quark or antiquark in that hadron. We first evaluate 
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the eikonal functions x(b, R(l, 1), R(2, 1)) in eq.(^) for the confining case, namely k — 1. 
The integration surfaces S\ and S 2 of the eikonal function 

Xp(M)^(2,l)) = (-^) 2 / dE^s) / dE"V) (F%(x, w)F%{x\ w)) A (49) 

are represented in figs. 3.3 and 4.1. The first of these figures gives a somewhat tilted three- 
dimensional view, while the second shows a projection on the transverse plane. The vectors 
Q(K,L) in the transverse plane connect the reference point C (with coordinates w) to the 
positions of the quarks and antiquarks of the loops 1 and 2. The quantity ip(K,L) is the 
angle between Q(1,K) and Q(2,L). 



In the integrations indicated in eq.([49"D we first note that the contributions involving front 
planes (Cyl 11 yl 12 , CBnBi 2 , CA 2 iA 2 2 and CB21B22) in fig. 3.3 vanish in the limit T — > 00, 
and therefore we are left with four remaining terms, the integrals over the products of the 
side planes of two different pyramids. The four- dimensional integration over the two surfaces 
can be finally reduced to a single integration^. The integrations along the directions x + and 
x_ can be performed, and result in expressions involving two-dimensional inverse Fourier 
transforms of the correlator D(k 2 ). A typical resulting contribution that comes from the 
product of surfaces (CAuBn) and (CA21B21) in fig. 3.3 is 

[da f 1 d(3 cos (1,1) {-\aQ(l,l) — /5Q(2, 1)| 2 ) , (50) 
Jo Jo 



where is the above mentioned two-dimensional Fourier transform of the correlator with 



n = 4 (see eq. (|37D ), defined as 

(-|f| 2 ) = _ 1 f d 2x QA^\K ± \ 2 g -j 

where £ is a two-dimensional vector of the transverse plane. This quantity, which is evaluated 
in Appendix 1, can be can be written in the convenient form 

^ 2 (4) (-|eT) = -^A 2 [(p|e|) 3 ir 3 (p|£|)] , (52) 

where 
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3tt 



(53) 



K3 is a modified Bessel function, and A2 is the Laplacian in two dimensions. 

Taking advantage of the Laplacian form, we can apply Gauss' theorem in two dimensions 
and eliminate one further integration^. 

It is useful to introduce a reduced eikonal function and a reduced loop-loop scattering 
amplitude through 
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K(g 2 FF) 



and 



JMM'(b, R\, R2) = 



WFF)\ 
(see eq.(^)). We then obtain 



Ri, R2) 



144 • 8 ■ iV^jV 2 - 1) 



(54) 



(55) 



X{b, R(l, 1), R(2, 1)) = - cos V(l, 1) I[Q(1, 1), Q(2, 1), ^(1, 1)] 

-cosV(2,2)/[Q(l,2),Q(2,2),^(2,2)] 
+ cos^(l,2)/[Q(l,l),Q(2,2),^(l,2)] 
+ cosV(2,l) J[Q(1,2),Q(2,1),V(2,1)] 



(56) 



where the quantities / are given by 

32 /3-7r^ 2 

I[Q(1,K),Q(2,L)M^L)) = —^ T ^ 

x{Q(l,AT) / [g(l,X) 2 + x 2 -2xg(l,iT)cos^(^,^)] 

JO 

■3tt 



— JQ(1, Kf + x 2 - 2xQ(l, K) cosip(K, L) 



dx 



+Q(2, L) / [Q(2, L) 2 + a; 2 - 2xQ(2, L) cos ^(iT, L)\ 
Jo 

■3tt 



i^ 



8 



Q(2, L) 2 + a; 2 - 2xQ(2, L) cos V (if, £) 



with Q(if, L) = \Q(K,L)\. 



(57) 
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From the eikonal function \ we contract the loop-loop amplitude J MM '(b, Ri, R 2 ) fol- 
lowing eq. (|3T)|) , where Ri and R 2 are shorthand notation for R(l, 1) and R(2, 1) respectively. 
The meson-meson scattering amplitude is then constructed by averaging over the transverse 
wave-functions, according to eq. (|4l|) . 

These results apply equally well to meson-baryon and baryon-baryon scattering if the 
baryon is represented as a meson-like structure in the quark-diquark picture. 

The evaluation of the eikonal function and of the observables for the non-confining cor- 
relator follows the same lines. Since this part of the correlator is not going to be used in 
the phenomenological analysis, we only present in the Appendix 2 the comparison between 
some corresponding quantities for the two cases, in order to exhibit their large differences. 

In order to treat the baryon as a genuine three-body configuration, we have to start from 
the expression given in eq. (|4"5|) . The projection of the loops for meson-baryon scattering 
is shown in fig. 4.2. In this case the relation between the scattering amplitude J and the 
eikonal functions x is more complicated than eq.([3~6"l). Let us now define x by 

A K = x(b, R(l, K), R(2, 1)) = - cos^(K, 1) I[Q(1, K), Q(2, 1), i/;(K, 1)] 

+ cos^,2) I[Q(1,K),Q{2,2),1>{K,2)] , (58) 

with the functions / given by eq. fl5"T|). The index K = 1,2,3 refers to the three quarks 
in the baryon loop, which is here taken as hadron number 1. As before, in the meson loop 
(hadron number 2) the quark is labelled 1, the antiquark is labelled 2. The relevant quantity 
for the transition probability is then 

J BM [b, R{1, 1), R{1, 2), R{1, 3), R{2, 1)] = 
' (A* + A\ + A\ - A 2 A 3 - A,A 3 - A,A 2 ) . (59) 



144 • 8 • Nq{Nq - 1 

For the scattering between two baryons in the three-body picture we define x by 

X(b, R(l, K), R(2, L)) = A KL = - cos^(K, L) I[Q(1, K),Q(2, L), i[)(K, L)] , (60) 

which obviously represents the contribution due to the interaction between the quark loop K 
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in baryon 1 and the quark loop L in baryon 2. Then the hadron-hadron scattering amplitude 
is given by 

J BB [b, R(l, 1), R(l, 2), R(l, 3), R(2, 1), R(2, 2), R(2, 3)] 1 



144 • 8 • iV^(JVg - 1) 
A\ x + + A 2 l3 + + A\ 2 + + + A 2 32 + A33 

— A12A13 — A22A23 — A32A33 — A11A13 — A21A23 — A31A33 

—A\\A\2 — A21A22 — ^431^32 — ^421^31 — A22A32 — ^23^33 
— A11A31 — Ai 2 A 32 — A i3 A 33 — AuA 2 \ — A12A22 — Ai 3 A 23 
+ ^ (^22^33 + ^23^32 + ^21^33 + A 23 A 31 + ^21^32 + ^22^31 

+A 12 A 33 + A 13 A 32 + A U A 33 + A 13 A 31 + A n A 32 + A 12 A 31 
+A 12 A 23 + A l3 A 22 + A n A 23 + A l3 A 2l + A U A 22 + A 12 A 21 



(61) 

In the simplest case we place the quarks around the central point of the baryon in such 
a way that the vectors R(l, K) form angles of 120 degrees (this configuration minimizes 
the string tension), and also choose the distances to the baryon center to be all equal, i.e. 
\R(1, 1)| = \R(1, 2)| = \R(1, 3)|. In this way, to form the amplitude, the configuration of the 
loops forming a baryon is specified by only one transverse vector, say #(1,1) of quark 1. 

For the hadron transverse wave-function we make the ansatz of the simple Gaussian form 

MR) = V^^-exp (-R 2 /S 2 H ) , (62) 

<->H 

where Sh is a parameter for the hadron size. 

Analogously to eq. fl41|) , we write the reduced hadron-hadron amplitude as an average 
over the hadronic wave-functions 

J Hl H 2 (b,SuS 2 ) = J d 2 R! J d 2 R 2 J Hl H 2 (b,RuR2) \MRi)\ 2 IM^f , (63) 

which is a dimensionless quantity. 

The hadron-hadron scattering amplitude in the eikonal approach is then given by 

T HlH2 = t2s[ K (g 2 FF)} 2 a 10 J d 2 b exp {iq- b) J Hl H 2 (b, S u S 2 ) , (64) 
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where the impact parameter vector b is in units of the correlation length a, and q is the 
momentum transfer projected on the transverse plane (in units of 1/a, so that the momentum 
transfer squared is t = — \q\ 2 /a 2 ). The eikonal approach requires large s and \t\ << s. 

We have verified that the contributions of next order in (g 2 FF) are smallil, which is a 
consequence of the presence of the colour factor in eq. (|36|) . 

Our normalization for Th x h 2 ls such that the total cross-section is obtained through the 
optical theorem by 

a T = ~ Im T Hi h 2 , (65) 



and the differential cross-section is given by 

da e£ 1 



T Hi h 2 \ 2 • (66) 



dt 167TS 2 

For short, from now on we write J(b) to represent J(b, S\, S2). 

The shapes of J(b) for the three cases of hadronic scattering are shown in fig. 4.3, 
against the impact parameter b (in units of the characteristic length a). In all the three 
curves represented in the figure, we have used S/a = 4 for both interacting hadrons. The 
label M means meson-like structure, while B means a three-body configuration. 

In all cases, J(b/a) can be written, in very good approximation, as a function of the form 

JiP) = J(0) %*+ff a f exp (-Ps(b/af) , (67) 

where Pi, P2, P3 are parameters, determined by fitting the exact (numerically obtained) 
values of J{b/a). 

Let us define the dimensionless quantities (as before, with b in units of the correlation 
length a) 



d 2 bb k J(b) , A: = 0, 1,2, ... (68) 
which depend only on S/ a, and the Fourier-Bessel transform 

I(t) = J d 2 b J {bay/\t\) Jib) , (69) 
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where J (bay\t\) is the zeroth-order Bessel function. Thus 

T HlH2 =is[K(g 2 FF)] 2 a 10 I(t) . 

Since J(b) is real, the total cross section a T , the differential elastic cross-section and the 
slope parameter(slope at t — 0) 



5 = ^fln 



dt V (it 



t=o 



(70) 



are given by 



/ = /, 6 tfFF)V, (71) 



and 



B = I h a 2 = Ka 2 . (73) 

We have here defined 

K = \ 'f , (74) 

which is a function of Si /a and S^/a only. 

We observe that in the lowest order of the correlator expansion used here, the slope 
parameter B does not depend on the value of the gluon condensate (g 2 FF). 

The curves for I = a T / [(Kg 2 (FF)) 2 a 10 ] and for K = B/a 2 as functions of S/a are 
shown respectively in figs. 4.4 and 4.5, for the BB, MM and MB cases, with hadrons of 
same extensions, S± — S2 — S. 

In the interesting ranges 1 < S/a < 3 and 0.5 < S 2 /S 1 < 1, parametrizations for the 
dependence of the total cross-section and slope parameter on the hadron size parameters Si, 
to an accuracy better than 3% are 

h=a {—) (75) 
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and 

B = 1.558a 2 + ( 7 /2) [S 2 + S 2 ] , (76) 
where the value of the constants a, (3, 7 are given in table 4.1 for the three different cases . 



Table 4.1 Values of the parameters to be used in eqs. (|75|) and ( [T6f) to determine total 
cross-sections and slopes. 



Case 


a 


P 


7 


MM 


0.00626 


3.090 


0.366 


BB 


0.00881 


3.277 


0.454 


BM 


0.00682 


3.135 


0.348 



A more precise parametrization, accurate in ranges of S\/ a extending from to 5, is 
given by 

I = a(-) P (77) 



and 

K = 1.558 + l(^) 5 , (78) 

where now the coefficients a and 7 depend on the ratio Si/ S\. The values of the parameters 
a, P, 7, S are given in table 4.2, for values of the ratios S2/S1 which we consider in the next 
section, where, besides pp and pp scattering, we also discuss pir, pK and pE scattering, in 
which cases Si 7^ S\. 



Table 4.2 Values of the parameters for eqs. ( [77]) and ([78|), for several values of Si/ Si 
which are important to represent the total cross-sections and slopes for different hadronic 
systems. 
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Hadrons 


1 


2 


parameters for cross- sect ions 


parameters for slopes 








a 


P 


7 


5 


5*2 = 5*1 


M 


M 


0.006260 


3.090 


0.3616 


2.023 




B 


M 


0.007846 


3.135 


0.4311 


1.955 




B 


B 


0.008814 


3.277 


0.4891 


1.892 


S 2 = 0.94,Si 


M 


M 


0.005610 


3.090 


0.3403 


2.023 




B 


B 


0.008000 


3.277 


0.4636 


1.892 


S 2 = 0.77,Si 


M 


M 


0.004159 


3.090 


0.2908 


2.023 




B 


M 


0.004446 


3.277 


0.3577 


1.955 


S 2 = 0.67,Si 


M 


M 


0.003353 


3.090 


0.2672 


2.023 




B 


M 


0.003604 


3.277 


0.3330 


1.955 



The t-dependence of the logarithmic slope of the differential cross-section is given by 

d A da eE \ 2 dl(t) 



B 



dtV" dt J I{t) dt 
a 2 r°° 27r 
t\ 



/ -=J{b)tfJr{baJ\t\)db 
Jo J\t\ 



I(t) 



For small t values, we obtain, expanding the Bessel function, 



(79) 



B(t) ~ B(0) 

L O 

= B(0)[1 + Ca 2 \t\] . 



l/27 2 h\ > I, 



(80) 



The values for C that we have obtained in our calculations for systems with S 2 — Si — S 
are almost independent of S/a. The values are C ~ —0.54 for MM systems (here is included 
the case of diquark picture for baryons) and C ~ —0.42 for BB systems (three-body picture 
for hadrons). 

5. CHOICE OF PARAMETERS AND COMPARISON WITH EXPERI- 
MENT 
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The numerical parametrizations of the total cross section and of the slope parameter B, 
eqs.(^) and (f76|) , are very convenient for comparison of our model with experiment. We 
first concentrate on elastic pp and pp scattering. In these channels data are available over 
a wide energy range. Since phenomenologically the Regge-pole parametrization works very 
welli'i, the vaccuum exchange contribution, to which our models refers, can be extracted as 
the pomeron contribution in a Regge-pole analysis. Donnachie and LandshoflS found that 
the parametrization 



rf om (PP,PP) = 21-70 mbs - 0808 (81) 

(with s in GeV) works very well over the whole range of data from yfs = 5 to 1800 GeV, so 
that we can use this expression for the pomeron contribution. According to our convention 
we choose as energy ^/s = 20 GeV, but we wish to enphasize already here that the value of the 
QCD parameters, namely the gluon condensate (g 2 FF) and the correlation length a entering 
the expressions ( [TT]) and (^) for the fundamental correlator are practically independent of 
the choice of the energy value. We return to this point at the end of this section. 

For the value of the logarithmic slope of the elastic cross-section at t = for -^i = 20 GeV 
we use@ 

B = 12.47 ± 0.10 GeV" 2 . (82) 

This value extrapolates very well, through a Regge amplitude with a'(0) = 0.25 GeV -2 , to 
the observedli value 5 = 17 GeV" 2 at v /i = 1800 GeV, so that it can be taken confidently 
as representative for the vacuum exchange contribution. 

Once the form of the correlator is fixed, we have two parameters in the model, which 
are fundamentally related to QCD, namely the gluon condensate (g 2 FF) and the correla- 
tion length a, and the extension parameter Sh, determined by phenomenological hadronic 
physics. Of course none of these parameters is completely free, as they appear in other phe- 
nomena besides soft high-energy scattering, and they may be obtained also through lattice 
calculations. In this section we take advantage of this independent information to check the 
consistency and fix the values of our physical parameters. 
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The gluon condensate < g 2 FF > was first determined in 1979 by Shifman, Vainshtein 
and ZakharovEl in the framework at QCD sum rules applied to the charmonium system, 
with a value (g 2 FF) = 0.47 GeV 4 . The analysis has been repeated and extended many 
times, several authors finding a range of values extending to considerably larger values (see 
Appendix 3). So we here accept as conservative estimate the range 

(g 2 FF) = (0.5 . . . 1.5) GeV 4 . (83) 

There are also theoretical uncertainties. Low energy theorems^ indicate that, in a world 
without light quarks, the condensate value would be larger by a factor 2 to 3 than the 
empirical value of our world with three light flavours. Therefore in our calculation we must 
refer to the value of the gluon condensate which is valid for a world without light quarks. 
Our model does not include dynamical effects of light quarks, which are treated as external 
sources moving on given (light-like) paths. A consequence of the absence of dynamical 
fermions is that the fermion (or Wilson loop) renormalization constant is equal to 1. Taking 
into account dynamical fermions not only would change the value of the gluon condensate 
but also would lead to a loop renormalization constant (see eq.QZID), 

" (k tlW{0 ' n) )A' (84) 



that is smaller than one. Thus in our case we should use a value for the gluon condensate 
from a pure gauge theory 

(g 2 FF) = (1 ... 3) GeV 4 . (85) 
As already mentioned in sec. 2, the fundamental correlator 

0) F£(0, 0)> A 



has been calculated in a pure SU(3) lattice gauge theory by Di Giacomo and Panagopoulosa 
using the cooling method. By this method the high frequency contributions are frozen 
out and therefore this correlator is just the one which can be compared to that of our 
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investigations. The cooling method works very well for Euclidean distances above about 
0.4 fm and the results obtained in the lattice calculations show that the confining tensor 



structure {8^5^ — 8^ a S vp } in eq.flTTD is definitely present and even dominant, with k ~ 3/4. 
In the physical range from r = 0.5 to r = 0.8 fm, the scalar function n(g 2 FF)D(—r 2 /a 2 ) is 
given by the function 

K(g 2 FF}D(-r 2 /a 2 ) = 24 Aexp(-r/X) , (86) 

with 

A = 3.6 10 8 A 4 , A = 1/(183A L ) and A L = (0.005 ± 0.0015) GeV . (87) 

Within the given accuracy, the scalar function D\ is proportional to D, and k ~ 3/4. 
Our choice of the correlator function D(z 2 /a 2 ) given by eq. (|37|) with n = 4, in the Euclidean 
region of the lattice results, namely for values 1 < —z 2 /a 2 < 9, is well approximated 
by an exponential function. We can therefore determine the gluon condensate and the 
correlation length a by fitting our expression for the correlator to the one obtained in the 
lattice calculation. The values obtained for k < g 2 FF > and a depend strongly on the value 
of the lattice-QCD parameter A^. In fig. 5.1 we show our correlator (solid line) together 
with the result found in the lattice calculation in the region 0.5 < r < 0.8 fm, for the choice 
A L = 0.0044 GeV. This fitting leads to the values k < g 2 FF >= 1.774 GeV 4 and a = 0.350 
fm. 

With this value for a, our correlator passes through zero for r m 1.4 fm. This change 
of sign is certainly an artefact of our special ansatz, but is has no practical consequences, 
since it occurs in a region where the exponential damping makes its contribution irrelevant 
anyhow. 

As mentioned in sec. 2.2, the evaluation of the (Euclidean) Wilson loop in the model of 
the stochastic vacuum yields a relation involving the condensate K(g 2 FF), the correlation 
length a and the string tension p (see eq. flTo]) ) . For our family of correlators the integral can 
be performed analytically (see eq . (|40|) ) and we obtain for the case n = 4 
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<9 2 FF) = 8 ^p . (88) 

Finally, we remark that the extension parameter S p is not completely arbitrary, but 
should be in a range of values S p = 0.5 ...1 fm, i.e. around the proton electromagnetic 
radius. 

We now use all this information to analyse our results . 

In fig. 5.2 we display the relation between the gluon condensate < g 2 FF > and the 
correlation length a obtained from different sources. Fig. 5.2a refers to the diquark picture 
and fig. 5.2b to the three-body picture for the proton. The solid lines show the possible 
choices of < g 2 FF > and a as obtained from our model, using eqs.([n]),([73]), (|75|) and ( |76|) 



and the experimental inputs <y pom = 34.9 mb and B = 12.47 GeV -2 . To indicate the ranges 
of values of the proton radius which are represented in the plots, we mark on these curves 
the points where the values of S p are 0.8 fm (3-body case), 0.9 fm (diquark case) and 0.6 fm 
(both cases). 

The dashed lines represents the results of the lattice calculation^, where the largest 
error comes from the uncertainty in the lattice QCD parameter = 5 ± 1.5 MeV. The 
values corresponding to some chosen values of are marked on this curve. The points for 
this curve have been obtained by fitting the lattice results to our form for the correlator, as 
exemplified for a given value of A^ in fig. 5.1. 

The dotted lines represent the relation between the gluon condensate, the correlation 
length and the string tension as obtained in the model of the stochastic vacuum00; for our 
form of correlator, this relation is given by eq . (p8|) . The upper and lower curves correspond 
respectively to string tension values p = 0.18 GeV 2 and 0.16 GeV 2 . 

As can be seen from the figures, the constraints from these three independent sources 
of information are simultaneously satisfied in a narrow region, providing a very consistent 
picture of soft high-energy pp and pp scattering, for the following sets of parameters. 

1) In the diquark picture 

a = 0.350 fm ; < g 2 FF >= 2.39 GeV 4 ; S p = 0.835 fm ; 
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A L = 4.4 MeV ; p = 0.18 GeV 2 . (89) 

2) In the three-body picture 

a = 0.361 fm ; < g 2 FF >= 2.08 GeV 4 ; S p = 0.730 fm ; 

A L = 4.2 MeV ; p = 0.16 GeV 2 . (90) 

The calculations described above lead to a determination of k < g 2 FF > . To obtain 
the value of the full condensate quoted above, we have used k = 0.74, as determined by the 
lattice calculation!!. 

The (pure gauge) gluon condensate is well compatible with the canonical value (see ^5|) . 
The resulting proton size parameter S p comes out quite close to the electromagnetic radiusEll 
value, which is R p = 0.862 ± 0.012 fm. The lattice parameter and the string tension p 
are also in their acceptable ranges. 

A very specific feature of the model of the stochastic vacuum is the dependence of the 
total cross-section on the hadron size, even if the latter is large as compared to the correlation 
length. The size dependence can best be tested by comparing the cross-sections for different 
hadronic systems. The Donnachie-Landshoff parametrizationl gives for the pomeron parts 
of pn and pK cross-sections the ratios 

<V/Opp = °' 63 ! VpK/Vpn = 0.87 . (91) 

In our treatment the theoretically predicted cross-section ratios depend on the hadron sizes. 
We take for the ratios of these sizes the ratios of the respective electromagnetic radii. The 
known values^ for these radii are 

R p = 0.862 ± 0.012 fm , R w = 0.66 ± 0.01 fm,R K = 0.58 ± 0.04 fm. (92) 

With this input we obtain for the predicted ratios the results of table 5.1 . 



Table 5.1 Ratios of the pomeron exchange contributions to the total cross-sections for 
different processes. 
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Cross-section 


p picture 


p picture 


Experimental 


ratios 


diquark 


3-body 


values 




0.66 ±0.02 


0.50 ±0.02 


0.63 




0.82 ±0.08 


0.82 ±0.08 


0.87 



We see that in the diquark-picture both ratios agree very well with experiment. In 
the three-body picture the <J p7T /<J pp ratio comes out too small. Presumably the assumption 
that the ratio of hadron radii is the same as the ratio of the electromagnetic radii is an 
oversimplification if the assumed structures for the mesons and baryons are so different, as 
when we relate a pn with a pp with protons in the three-body picture. This is confirmed by 
the fact that both pictures reproduce perfectly well the a pn /a P K ratio, where the effect of 
the three-body baryon structure is the same in both terms of the ratio, cancelling out. 

With the criterium for the ratios of hadronic radii now fixed, we can also calculate the 
7i7i cross-section and we find that for both the diquark and the three-body picture the 
factorization relation 

is numerically perfectly fulfilled, as can be verified through the parameter values given in 
table 4.2. From the data on the pp and p7i systems we thus predict the pomeron part of the 
7i7i scattering cross-section to be 

= 8.6 mb s - 0808 . (94) 

The only experimental result, known to us, on the tttt cross-section@ gives cr W7r = 10 mb 
at y/s = 4 GeV, while the above expression predicts 9.6 mb for the pomeron part at this 
energy. 

There are some data available on hyperon-proton scattering)!, which also fit nicely in the 
general picture. If we dare to extract a vacuum exchange part from the few existing data 
for Ep scattering, we obtain 
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a pS - = 19.6 mb s 0808 . (95) 

From this cross-section we obtain for the ratio of the extension parameters, in both the 
diquark and three-body pictures, the value 

S^/S p = 0M, (96) 

which certainly is a reasonable result. 

The size dependence of the slope parameter at t = is given in parametrized form 
by eq. ([76]) . In table 5.2 we show the results of our model for the differences of the slope 
parameters for pp , T,p , pir and pK scattering. The extension parameters Sh for the different 
hadrons are assumed to be proportional to the electromagnetic radii given in eq . (|9~2|) . The 
experimental numbers for the vacuum exchange part of the slopes were taken from the 
analysis by Burqi (see their table 7, extrapolated to t — according to table 8). The 
extrapolation to t — has little effect on the value of the slope, for our purposes at this 
point. 



Table 5.2 Differences of the logarithmic slopes of the pomeron part at t = for different 
processes, in GeV~ 2 . 



Slopes 


Diquark 


3-Body 


Experimental 




picture 


picture 


values 


Bpp — B^p 


0.40 


0.45 




Bpp B-xp 


1.30 


1.45 


2.48 


B n p — B Kp 


0.43 


0.34 


0.34 



Comparison with the data shows that the difference of the slopes for irp and pp scattering 
is underestimated in our model, while the difference for np and Kp is well compatible with 
experiment. 
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The parametrization (|7q) of the slope parameter 

B = 1.858 a 2 + 0.183 (S 2 + Sf) , (97) 

where Si and S% represent the hadron sizes, and a is the correlation length, is a good 
approximation, within a few percent, to the results of our model. This parametrization can 
be compared to that of a modified Chou-Yang pictureil, which allows also for a quark form 
factor 

B = ^(R 2 q + R\ + Rl) , (98) 

where R q and Ri are the electromagnetic radii of the quark and hadrons respectively. Since 
in our results 1.858a 2 is about 6 GeV -2 , we see that in our model the correlation length a 
gives a much more important contribution to the slope than the quark-form factor of the 
modified Chou-Yang model. These different predictions should be tested experimentally. 

It is a curious result of our treatment that, given the ratios of the electromagnetic radii of 
the hadrons, the ratios of the differences of slope parameters are practically those of simple 
integers. Thus in the diquark picture 

9 

B p7T — B pK : B pp — B p7T : B pp — B pK : B pp — £> pS = 1:3:4: — - (99) 

We remark that the last figure in this sequence of ratios is a bit more uncertain than the three 
first ones, because it was determined using the poor scattering data, and the predicted, 
not experimentally measured, ratio (S^/S p ) = 0.94. 

Our model also predicts a t-dependence for the logarithmic slope parameter, as shown 
in eq. flSC| . In the analysis by Burqi, the t-depencence of the elastic cross-sections has been 
written in the form 



da da 
dt dt 



e [B(0)+ct+dt 2 ]t ^ ^qq^ 

t=0 



with the parametrization for B(t) 

B{t) = B(0) - 2c\t\ + 3d\t\ 2 . (101) 
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From eq. ([80|) and our results for the slope of the slope C, we obtain for this parameter c the 
values 

c = 10.6 GeV -4 in the diquark picture , (102) 
c = 8.8 GeV -4 in the three-body picture . (103) 

These values are to be compared with the experimental value 6.8 ± 0.5 GeV" 4 obtained by 
BurcS. However we should not overestimate the significance of this comparison. On the 
theoretical side this slope of the slope depends strongly on the precise form of the profile 
function J(b), while the experimental results are not accurate and may still be contaminated 
by Coulomb interference. 

If we consider hadron-hadron scattering for two hadrons of equal sizes, we can eliminate 
the hadron radius between the parametrized forms for a T and B, and obtain the relation 

T = ^-P/Ss . Jpp . X2„ 10-2/3/5/ R _ , o R o„2\/9/* 



(T P om = a l ' 



< g 2 FF >) 2 a 10 - 2 ^\B - 1.858a 2 ) 13 / 5 , (104) 



where the values for the model dependent parameters a, (3, 7, 5 are given in table 4.2 . 

If we fix n(g 2 FF) and a at one given energy (e.g. ^fs = 20 GeV), we obtain through 
eq. (|104|) a parameter- free relation between the total cross-section <rJ om and the slope param- 
eter B. 



In fig. 5.3 we display cxJ om against B as given by eq. (|104|) , for both diquark and three- 
body pictures for the proton, using the sets of parameters given in table 4.2 . In the same 
figure we also show the relation 

4e S9 e = ^ e^*-^ , (105) 

obtained from a Regge amplitude using the slope of the pomeron trajectory a'(0) pom = 
0.25 GeV" 2 , and as input at y/s = 20GeV the values = 35 mb, and B = 12.47 GeV 2 . 
The experimental data! at 540 GeV and 1800 GeV are marked in the plot, together with 
the input data at 20 GeV. We observe that our relation ( |104| ), which contains no free 
parameters, describes the experimentally observed relation between B and a T astonishingly 

36 



well. Besides, there is a surprising agreement between the Regge parametrization line and 
our results. It must be remarked that the constant term 1.858 a 2 in our expression for the 
slope B is important for this good agreement with experiment. 

The application of our results to different energies implies a very slow energy dependence 
of the hadronic radii. An explicit relation is obtained if we bring into eqs.(|7ID and ( [77| ) a 
parametrization for the energy dependence of the total cross-sections, such as the Donnachie- 
Landshoflfl form of eq.(fy]). We thus obtain for the proton radius 

^0.0808/3.090 



S p (s) 



K(g 2 FF} 2 a 10 (21.7mh)i 1 '' , " ,il 



0.00626 



in the diquark picture, and 



S p (s) 



K(^FF) 2 a 10 (21.7mb)lV3.277 ^ ao8o8/ ^ 77 



s U.U8U8/3.277 (107) 



0.00881 

in the 3-body configuration case. The values thus obtained for S p are in the region of the 
proton electromagnetic radius. 
6. CONCLUSIONS 

We have obtained a very consistent description of the data on soft high- energy scattering. 
Our basic assumption is that the low frequency, i.e. nonperturbative, contributions to the 
scattering amplitudes can be approximated by a Gaussian process. The parameters deter- 
mining the observable quantities are the gluon condensate and the correlation length of the 
vacuum field fluctuations. These parameters occur also in completely different connections 
(SVZ sum rules, lattice calculations, low energy hadron spectroscopy), and all the conditions 
posed on them can be consistently satisfied. The third parameter entering our calculations 
is the transverse hadron size, which may be related to the electromagnetic radius. 

Gauge invariance is observed on each step of the calculation, which is based on a formal- 
ism for loop-loop scattering, rather than on a quark-quark scattering picture. In this way 
even finite distance correlations of the vacuum field tensor lead to long range correlations 
that are the common source of confinement, and to a dependence of the total scattering 
cross-section on the hadron size. This mechanism can be interpreted as a string-string 
interaction in hadronic scattering. 
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The size dependence of the cross-sections leads to a natural explanation for the exper- 
imentally observed flavour dependence of the total cross-section. The ratio between the 
pion-proton and (anti)proton-proton cross-sections emerges in our model as a consequence 
of the different hadron sizes, and the factor 2/3 comes out as a consequence of the ratio of 
the electromagnetic radii, and not from quark additivity. 

The size dependence of the slope parameter, i.e. the logarithmic slope of the differential 
elastic cross-section, can be parametrized in a form similar to that of a modified Chou-Yang 
rnodel^ with a finite quark radius, which in our model appears as a correlation length. 

Elimination of the extension parameter yields a parameter-free relation between the total 
cross-section and the slope of the elastic cross-section which agrees very well with experiment. 

The investigations described in this paper can be extended in many directions. In the 
present calculations only one of the two possible tensor structures determining the low 
frequency contributions is taken into account. The inclusion of the second term, which 
could also describe perturbative effects (and even Coulomb interaction), would pose no 
important technical problems. Furthermore, we have restricted ourselves to the lowest order 
non-vanishing contribution, which is quadratic in the gluonic correlator. We have checked 
that this is justified for the total cross-section and the slope parameter. But our amplitude is 
purely imaginary, and quantities like the p-parameter (the ratio of the real to the imaginary 
parts of the elastic scattering amplitude) can only be described if we go one further order 
in the contributions to the correlator. 

A further important step to be developed is the test of the importance of the factorization 
implied by the assumption of a Gaussian process. This is certainly crucial in the present 
investigation, with possible consequences for the phenomenological analysis, but it remains 
to be seen which of the more general features depend on this approximation. 

It would of course be highly desirable improve the present model in order to describe 
the dependence of the cross-section with the energy. However, although we cannot obtain 
the (slow) rise of the total cross-section with the energy, we can eliminate the hadron radius 
parameter S and obtain a parameter free relation between the total cross- section and the 
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slope parameter B, as shown in fig. 5.3 . 
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Appendix 1 . THE CORRELATION FUNCTIONS 

In eqs.flllD and fl28|) there appear two independent arbitrary scalar functions, D(z 2 /a 2 ) 
and Di(z 2 /a 2 ), which are supposed to fall off at large distances with characteristic lengths 
a, called correlation lengths, and must have forms that can be analytically continued from 
Euclidean to Minkowski space-time descriptions of field theory. If the expression for the 
correlator is used in Euclidean QCD, the scalar function D(z 2 /a 2 ), which is zero in Abelian 
theories (if there are no monopoles) leads to a linearly rising potential, namely to confine- 
ment. QCD lattice calculations have shown that the dominant contribution to the correlator 
actually comes from the term with D(z 2 /a 2 ) , namely k ~ 1. Besides that, in Appendix 
2 we show that the non-confining part Di of the correlator has much less influence on the 
values of the eikonal functions of high-energy scattering. Consequently, in the present work 
we take into account only the confining term D(z 2 /a 2 ), with the weight k determined by 
the lattice calculations, neglecting the effect of D\[z 2 /a 2 ) altogether. 

We thus concentrate on D(z 2 /a 2 ), and take as an ansatz the family of functions 

D<n)ie> = - m I Wf> (i^w exp Hfc?/A " ) ' " - 4 • (108) 

where 

£ = z/a, (109) 

a is the characteristic correlation length, and the constants A n and A n are to be fixed by 
normalization. It is convenient to absorb A n into k through k/X n — > k ; then £,k and A n are 
all dimensionless. In the Euclidean metric 

_ id A k = d A K , K 4 = ik , k 2 = —K 2 = -(\K\ 2 + K\) , (110) 

for space-like vectors 

m^^es), £ 2 = -!£?, (in) 

we have 
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£>(«)(£ 2 ) = DW(-|eT) = Q(-l) n+1 J 



d 4 K A n X e n K 2 ~ 

exp («A ■ 4) 



(2tt) 4 (A 2 ^ 2 + ! 

The constants A n and A n can be fixed by the normalization conditions 



;ii2) 



dW(o) = i , / dm rt n \-\tf) 

Jo 



(113) 



The second of these two relations has the role of a definition for the correlation length. 
These calculations can be made analytically, leading to 

4tt 2 



An = ("I) 



71+1 



and 



A, 



-(n-l)(n-2)(n-3) 



4 r(n-3) 



(114) 



p n 3v^r(n-5/2)- 
For simplicity of notation, from now on we use more often p n = 1/A n , instead of X n . 
The integrations in eq. (|112|) can be performed analytically. We obtain 



(115) 



£> (n) (-ieT 



■l) n+1 3A n -> 



n—3 



[n - l)K n ^{p n \£\) - ±{p n \£\)K n _ 2 {p n \£\) 



(116) 



K 2 2 n - l T(n) 

which is the general form for the class of functions considered. K v (x) is the modified Bessel 
function. 

Alternative useful forms to eq.( |116| ) are 



B<«>(-|f 



ir 2 2 n T(n) 



d( P n\t\) 



(Pn\Z\V +1 K n _ 3 (p n \Z\) 



117) 



and 



7c 2 2 n T n) 



d(p n \t\) 



(Pn\Z\f 



d 



d(Pn\ti\) 



[(Pn\i\) n - 2 K n . 2 ( Pn \i\)] 



(118) 



These correlation functions are negative at large distances, behaving like 



2 n - 2 Y{n - 3) 



(Pn|el) n - 5/2 exp(-p n |e|) 



(119) 



After the limits are taken, which make the long sides of the rectangular Wilson loops 
tend to ±oo in the x% direction, the remaining variables in the integrands are coordinates 
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of points in the transverse plane. The distances z between such points enter in the final 
expressions for the eikonal functions % as arguments of the two-dimensional inverse Fourier 
transform of D(k 2 ) in eq. fl3~T|) , which is defined by 

7 (n)/ I 



6A n X 6 n \K ± \ 2 



(27r) 2 J - (_ A 2|^ ± | 2 -1)« 

where £ is any two-dimensional vector of the transverse plane 
Using the values for A n and p n written above, we obtain 



exp (iK ± ■ f ) 



(120) 



■{CI' 



2( 8 -»)r(n - 3) 



9[r(n-5/2)] 2 
Important alternative forms are 



(Pn^l) 



n-2 



[n - l)K n _ 2 (p n |£|) - -{p n \i\)K n ^{p n \i\) 



;m) 



in). 



2 (7-n)r( n - 3) 



9[r(n-5/2)] 2 
and 

7 {n)< 2( 7 -™)l>-3) 



(Pnlel)- 1 



d 



d(f>n\t\) 



( Pn iei)^- 2 (p n iei) 



(122) 



•TTX-ieP 



9[r(n - 5/2)] 



r(Pniei)" 1 



(PnlCl) 



d(j»n\t\) 



(p n \^\r-"K n ^(p n \^\) 



(123) 



or 



where 



(n). 



■leT 



2 (7-n)p( n _ 3) 

'9(r(n-5/2)) 2 



A 2 ^(P«lf1) , 



(124) 



V> (n) (Pn|£'|) = (^leD^^n-rCPnlel) , 



(125) 



and A 2 is the 2-dimensional Laplacian operator. This Laplacian form is important in our 
calculation, as it allows lowering the order of the integrations, through Gauss theorem. 

The moments of the D^(— |<fj 2 ) and T^^— \£\ 2 ) functions can be readily obtained from 
eqs.( |117| ) and ( |122| ) respectively. For D^ n > we obtain 



M p [D n {-\H a )] = I (Pn\Z\Y D( n \-md(p n \Z\) 

42 



(126) 



where we observe that A4s[D n (— 1£*| 2 )] = and that all moments higher than p = 3 are 
negative, for every n. Thus all correlation functions have a zero. 
For we have the moments 

M p \r?\-K\ 2 )) = 1°°(pM) p ■ffVlfTMft.kl) 
2(P+4 ':i - P ) p ''t r( "' s) r(e±±)r(n + !^*) . (127) 



9 v [r(n)] 2 V 2 J V 2 

Thus for every n the moment p = 1 vanishes, and all higher moments are negative. Thus 
all functions pass through zero. 

We must choose the value of n in order to fix a specific form of the ansatz for the 
correlation function. Actually, the dependence of the final results on this particular choice 
has been testedfH, and found to be not very marked. The reason is that all correlation 
functions are normalized to 1 at the origin, and decrease exponentially at large distances. 
It is enough that the chosen function falls monotonicaly and smoothly in the range of 
physical influence (up to about one fermi, say), and there cannot be much difference in the 
results obtained using different analytical forms. Of course there will be differences in the 
specific values given to the correlation length parameter a, due to the different values of the 
multiplicative factor p n in the argument of the exponential behaviour (see eq.( |119|) ), but 
such differences can be taken into account and absorbed when different forms of correlation 
functions are compared. 

In the present work we make the choice that n=4, which in the Euclidean region leads to 
a good representation of the lattice calculations!!. We then have for the correlation function 



£ (4) H£?) = NCI) 

and for the 2-dimensional inverse Fourier transform 



128) 



^ 4) (-ie?) = ||(p 4 iei) 2 



2Ao(/;,K|)- {-^-P4\Z\)l<df»\L\) 

P4Q 
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-^A 2 [(p 4 |e|) 3 ^ 3 (P4|el)] , (129) 



where 
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3tt 



(130) 



The functions L>( 4 )(-|£| 2 ), ^Hli 2 ) and ^ (4) (p 4 |^) are represented in figs. (7.1), (7.2) 
and (7.3), against the variable x = p^\. The correlation function D^\— |£| 2 ) has a zero 
at Pa\^\ = 4.43 while J r ^\— \^\ 2 ) has a zero at p^\ = 3.05. As we will see from our final 
results, the locations of these zeros are beyond the range of physical influence. 

Practical representations for these functions, that are important for the numerical work, 
can be obtained. As a tool to obtain the parameters in approximate representations, we 
may use the moments of the functions, which are explicitly given by eqs. (|126|) and (|127|) . 
The function 



^\x) = x 3 K 2 {x) 



;i3i) 



where x = /?4|£| , enters in heavy numerical computations to produce the eikonal functions. 
It appears in eq.(^) as 

-x 3 K 2 (x) . 



dip 
dx 



Let us write 



■*// 4) (x) ~ 8 [1 + x + a x x 2 + a 2 x 3 + a 3 x 4 ] e x , 
that satisfies the constraints at the origin ip(0) = 8 and ip'{Q) = . We then have 

2 8 r 



(132) 



9tt 



-2 + 4ai + 1 -5ai + 9a 2 )x 



+ (ai — 7a 2 + 1603)0; + (a 2 — 9a 2 )x + a^x 



Using properties of F^i— |£j 2 ), we obtain 



F?\x = U) = 2 ! /{^) 



— T 

dx 2 



(4) 



= 0) = 

r T^dx = 16/3 
Jo 

x dx = 



(ii = 3/8 , 

a 2 = 1/24 , 

17 3tt, 
a3 "8 ( 3~T' 

identity . 



(133) 



(134) 
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There are no more free parameters. The next moment 



Jo 



oo 



x 2 Tf ] dx = -80/3 



is reproduced in good approximation. 



The representations (|132|) and ( |133| ) are excellent for our purposes, and can be safely 



used, reducing substantially the computation time. 

Appendix 2. Influence of the confining and non-confining correlators on 
high-energy scattering 

We have seen in subsec. 2.2 that only one of the tensor structures present in eq.(^8|) 
leads to confinement, namely the part of the correlator whose scalar function is denoted by 
D. If k = 1 this is the only existing contribution, and since it has the property of leading 
to confinement, we refer to it as the confining case . The other tensor structure, which for 
k = is the only one contributing, is referred to as the non-confining case. Comparing the 
effects of the two functions, we put 



in order to have the same spatial behaviour for their contributions to the correlator. 



X is plotted as a function of b/a, i.e. the impact parameter in units of the correlation length 
a. In the confining case the eikonal function is approximately proportional to the overlap 
region of the two loops in the transverse plane. In the non-confining case there is a large 
contribution if both quarks or antiquarks from loop 1 and loop 2 coincide, i.e. for 6 = 0; the 
contribution is smaller if there is coincidence of only one particle (quark or antiquark) from 
of each loop. Thus in the confining case we have typically a true string-string interaction, 
while in the non-confining case we have a quark-quark interaction. 

In fig. 8.2 we show the strong dependence of the eikonal function \ on the orientation 
of the loops for the confining case. The figure shows \ as a function of b/a, for the case 




D ± (k 2 ) = D{k 2 ) . 



(135) 



In fig. 8.1 we represent the reduced eikonal function x(6, R(l, 1), R(2, 1)) for the case 
where b, R(l, 1), and R(2, 1) are all parallel and with a ratio \R\/a = 2. The eikonal function 
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that the loops are parallel in the transverse plane i.e. 9\ = 62 = 0. The common angle 6 is 
a parameter, running from 0° to 80° in steps of 20°. In order to strengthen and make very 
visible the effect, we have chosen a large value \R\/a = 10. From the figure we observe that 
the loops interact as one- dimensional objects in the transverse plane. 

In fig. 8.3 we show the cross-section (in arbitrary units) as a function of the hadron 
extension S/a, using a Gaussian wave function ^2/tx{1/ S) exp(— R 2 / S 2 ). For S/a < 1, the 
two cases nearly coincide, but for larger values the cross- section in the non-confining case 
stays nearly constant, whereas the cross-section for the confining case continues to increase 
approximately like (S/a) 3 . Thus we have quark additivity only for the non-confining case. 
Since the correlator for the confining case can only be present in a non-Abelian theory, an 
Abelian model will always yield quark additivity if (S/ a) > 1 holds§. 



Appendix 3. The value of the gluon condensate < g 2 FF >. 

The most accurate determination of the gluon condensate < g 2 FF > comes from sum 
rule analyses of the charmonium system. In their original paper on sum rules, Shifman, 
Vainshtein and ZakharovII obtained the value 

(g 2 FF) = 0.47 GeV 4 . (136) 

This analysis was extendecfi, yielding, taken into account all observed ground states of 
given quantum numbers in the charmonium system the range of values 

0.51 GeV 4 < (g 2 FF) < 0.79 GeV 4 . (137) 

Other analyseJ^H^ yield considerable higher bounds for the gluon condensate. The main 
uncertainties areS the error in the pole mass of the charmed quark, radiative corrections 
and the contributions of higher condensates. The last point is the most difficult to controlil. 
Models indicate that they could increase considerably the value of the gluon condensate0'@. 
A conservative estimate for the gluon condensate is 
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(g FF) = (0.95 ± 0.45) GeV 4 . (138) 

NovikovS has however argued convincingly against the larger values (say larger than 0.8 
GeV 4 ) for the reason that the analyses yielding large upper bounds for the gluon condensate 
either rely on only one specific channel, or are based on particular models. A non-diagonal 
sum rule analysis of the matrix element (OlA^lvr), with 

A M = \-g e PfMa(3 d 7 P F a/3 u , (139) 

that is particularly sensitive to the gluon condensate, yields also a low value. 

There is however another theoretical difficulty In a world without light quarks the 
value of the gluon condensate could be considerably larger. A low-energy theorem has been 
derived^ relating the change of the gluon condensate with respect to (light) quark masses 
with the quark condensate 

d ^ 967r 2 

< g FF >= — — <qq> , (140) 



dm q Pi 

where < qq > is the light quark condensate and f3i is the first coefficient of the Gell-Mann 
Low function. For SU(3) the factor is large, 96tc 2 / ' (3i ~ 100. 

Since the quark condensate < qq > is negative, an increase of the quark mass leads to 
an increase of the gluon condensate. In a world of pure gauge fields the gluon condensate 
would thus have a larger value. This pure gauge value is obtained from the empiric one, 
which refers to a world with three light flavours, by taking all masses to infinity, i.e. by 
integrating eq.( |140| ) up to a value of m q where the quarks decouple. The authors of the 
theorem^ estimate that the pure gauge value is about two to three times higher than the 
empirical value. As argued in sec. 5, in our model we must use the pure gauge value of the 
gluon condensate. 
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Figure Captions 

Figure 3.1 Wilson loops formed by the paths of quarks and antiquarks inside two mesons. 
The impact vector b is the distance vector between the middle lines of the two loops. R\ 
and i?2 are the vectors in the transverse plane from the middle lines to the quark lines of 
meson 1 and 2 respectively. For the antiquarks the corresponding vectors are —R% and — i?2- 
The front lines of the loops guarantee that the mesons behave as singlets under local gauge 
t r ansf or mat ions . 

Figure 3.2 Representation of the non-Abelian Stokes theorem. The contour integral 
running from w to im then to x^) an d back to w is deformed in order to become a surface 
integral. Here and x@) represent the coordinates of the quark and antiquark in a meson, 
respectively. 

Figure 3.3 Tilted perspective view of the surfaces S% and S 2 obtained from the line 
integrals along the Wilson loops dS\ and dS2 after appying the non-Abelian Stokes theorem. 
The line BnAn is the quark path line and i?j2^4«2 is the antiquark path line of meson i. C 
is the reference point (with coordinates w in the non-Abelian Stokes theorem). 

Figure 3.4 Wilson loops (without traces) describing a baryon path in the three-body 
picture. The line from to x',^ represents the path for the quark labelled i. The line from 
Xh to x' h is the path for the central point of the baryon. 

Figure 4.1 Geometrical variables of the transverse plane, which enter in the calculation 
of the eikonal function for meson-meson scattering. The points C\ and C2 are the meson 
centres. In the integration, P2 runs along the vector Q(2, 1), changing the length z, which 
is the argument of the correlator characteristic function. In analogous terms, points P\, Pi 
and P2 run along Q(l, 1), Q(l, 2) and Q(2, 2). This explains the four terms that appear in 
eq-© • 
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Figure 4.2 Geometrical variables for the calculation of the eikonal function for meson- 
baryon scattering. Notation analogous to fig. 4.1 . 

Figure 4.3 Dimensionless function J(b/a), where b is the impact parameter, and a is the 
correlation length. The three represented cases refer to baryon-baryon {BB), meson-baryon 
(MB) and meson-meson (MM) amplitudes, with the same values (S± = S% = S = 4a) for 
the hadron extension parameters. 

Figure 4.4 Function Iq(S/ a) representing the dependence of the total cross-section on the 
extension parameter S (in units of the correlation length a) of the hadron wave-function. M 
represents a meson-like hadronic configuration ( qq for mesons, quark-diquark for baryons), 
and B represents a 3-body star-like picture for baryons. We have here used Si = S 2 = S in 
all cases. 

Figure 4.5 Same as for previous figure, for the quantity K(S/a), which represents the 
slope parameter B divided by a 2 . 

Figure 5.1 Demonstration of the determination of the QCD parameters through the 
fitting of our correlator to the lattice calculation resultsil, for a given value of A^. Dashed 
line: lattice results for = 4.4 MeV. The arrows indicate the range inside which the lattice 
calculations were made. Solid line: our fitted correlator, with the fitted values a = 0.350 fm 
and K(g 2 FF)D(-\£\ 2 ) = 1.772 GeV 4 at = 0. 

Figure 5.2a Constraints on the values of K,(g 2 FF) and of the correlation length a. Solid 
line: Our model of high energy scattering, with a^ p = 35 mb and B = 12.47 GeV -2 , using a 
diquark picture for the proton ; dashed line : fit of our correlator to the lattice calculation^!; 
dotted lines : relation obtained from string tension p, eq. (p8|) , the upper curve corresponding 
to p = 0.18 GeV 2 , and the lower one to p = 0.16 GeV 2 . 

Figure 5.2b Same as fig. 5.2a, for the case of a 3-body picture for the proton. 
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Figure 5.3 Relation between the total cross section and the slope parameter B. Solid 



lines: predictions given by eq. ( |104j ), obtained from our model by eliminating the hadronic 
size, in the cases of diquark and 3-body pictures for the proton ; dashed line: relation 
obtained from the Regge amplitude with a(0) = 1.0808 and a'(0) = 0.25 GeV -2 ; star : 
input data at a/s = 19 GeV ; crosses : experimental data at 541 and 1800 GeV. 

Figure 7.1 Correlation function D^'(— |£| 2 ) , given by eq.( |128j ), against the variable 
x = p4^| = (37r/8) (z/a) , where z is the physical distance, and a is the correlation length. 
The correlation function is normalized to one at the origin. 



Figure 7.2 Two-Dimensional Fourier transform F^i - 1£*| 2 ) j given by eq.( |129|) , against 
x = p 4 |H = (3vr/8) (z/a) . 

Figure 7.3 Function ip^>(x) = x 3 K 3 (x) , given by eq.( |130| ), and represented approxi- 



mately through eq. ( |132| ) 



Figure 8.1 Comparison of the values of the reduced eikonal functions \ for the extreme 
pure confining (k = 1) and pure non-confining [k = 0) cases, in a situation where the three 
vectors b , R(l, 1) , and R(2, 1) are parallel. 

Figure 8.2 Angular dependence of the reduced eikonal function x for the purely confining 
case (k = 1). The figure shows results for configurations where the two interacting loops are 
parallel (6*i = 62 = 9) and \R\\ = I-R2I = 10a . The angle 9 varies from 0° to 80° in steps of 
20°. 

Figure 8.3 Comparison of the total cross-sections obtained in the two extreme conditions 
of K — 1 (pure confining case) and k = (pure non-confining case). 
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